On the Origin of the 2DEG Carrier Density at the LaAlOa/SrTiOa Interface 
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Transport measurements of the two-dimensional electron gas (2DEG) at the LaAIOa/SrTiOa in- 
terface have found a density of carriers much lower than expected from the "polar catastrophe" 
arguments. From a detail density-functional study, we suggest how this discrepancy may be rec- 
onciled. We find that electrons occupy multiple subbands at the interface leading to a rich array 
of transport properties. Some electrons are confined to a single interfacial layer and susceptible 
to localization, while others with small masses and extended over several layers are expected to 
contribute to transport. 

PACS numbers: 73.20.-r,73.21.-b 



In a seminal paper, Ohtomo and Hwang^ reported the 
observation of a conducting layer of electrons at the inter- 
face between two insulators LaAlOs (LAO) and SrTiOs 
(STO). The carrier density was quite high, of the or- 
der of 10^^ cm^^, as also the mobility, leading to the 
prospect for high-performance oxide electronics. Conse- 
quently, considerable attention has been focused on the 
understanding of the origin of these charge carriers/^ Sys- 
tematic experiments on oxygen annealed samples as well 
as on samples grown under high oxygen pressure con- 
ditions have tuned down the carrier density to about 
1 — 2 X 10^^ cm~^, which is now the common carrier den- 
sity observed by several groups i^'^'^'i Thus the emerg- 
ing consensus is that the higher carrier density ~ 10^^ 
cm~^, originally observed, is due to the electrons released 
by the oxygen vacancies and they have a 3D character, 
while the lower carrier density 10^^ cm~^ may be an 
intrinsic feature of the interface. However, this intrinsic 
carrier density is much smaller than the 0.5 electrons per 
unit cell (3.5 x 10^'* cm~^) needed to suppress the "po- 
lar catastrophe,"^ which arises because of the diverging 
Coulomb field due to the polarity of the individual lay- 
ers (Fig. [1]). Thus the origin of the carrier density for 
the intrinsic interface is poorly understood and remains 
a key puzzle for these interfaces. 

In this paper, we suggest a solution for this puzzle 
based on the results of our detail density- functional study 
of the electronic structure. Although there also exists a 
second type of interface, viz., the p-type interface with 
SrO/A102 termination, we focus on the n-type interface 
in this paper which is illustrated in Fig. 1. Our cal- 
culations reveal the presence of different types of bands 
at the interface. We argue that some of these electrons 
may become localized at the interface due to disorder 
or electron-phonon coupling and may not participate in 
the conduction, while the others form dclocalized states 
contributing to transport. 

Density functional calculations were performed us- 
ing the generalized gradient approximation (GGA) for 
the exchange and correlation.^ We have not included 
the Coulomb corrections within the LDA-I-U functional, 
as it may not describe the correlation effects prop- 
erly for low carrier density situation, away from half 
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FIG. 1: Nominal layer charges at the n-type interface with 
LaO/Ti02 termination. Bottom figure shows the diverging 
Coulomb potential ("polar catastrophe", solid line) and an 
example of how it may be avoided by adding half an electron 
to the single Ti02 layer at the interface (dashed line). 



filling. The Kohn-Sham equations were solved by ei- 
ther the linear augmented plane waves (LAPW)^^ or 
the linear muffin-tin orbitals (LMTO) methodsii. Ear- 
lier density-functional studies have been reported in the 
literaturei^ii^iiiii^; however, since the supercells used in 
these studies were quite small, the nature of the interface 
states could not be determined, which is the primary fo- 
cus of our paper to develop our arguments. 

We used the supercell (STO)7.5/(LAO)7.5 for studying 
the interface electron states, where the half layers in the 
formula signify an extra layer of Ti02 or LaO, so that we 
have two identical n-type interfaces (Ti02/LaO) present 
in the structure. The atomic relaxation was performed 
using the LAPW method for a slightly different super- 
cell, viz., for (STO)8.5/(LAO)4.5 to keep the calculations 
manageable. This is reasonable since results presented in 
Fig. [2] show very little relaxation in the LAO part be- 
yond the second layer, which in turn is consistent with 
the small charge leakage to the LAO side. 

The GGA optimized bulk lattice constants are 3.94 A 
and 3.79 A for the STO and the LAO bulk, respectively, 
which are quite close to their experimental values (3.905 
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FIG. 2: Computed relaxed atomic positions, consistent with 
the spread of the interfacial electrons into many Ti02 layers. 

and 3.789 A). For the optimization, we fixed the lattice 
constant for the Sr part to be 3.91 A, both in the direc- 
tions parallel and normal to the interface as appropriate 
for samples grown on the SrTiOs substrate, while the lat- 
tice constant for LAO normal to the interface was chosen 
to preserve the volume of the bulk unit cell. This defined 
the supercell which was held fixed, while the internal 
positions of the atoms were allowed to change in a se- 
ries of self-consistent calculations until the forces became 
smaller than about 6 mRy/ au. The main effect of the re- 
laxation, shown in Fig. [2l is to polarize the cation and the 
anion planes near the interface similar to the relaxation 
for the well-studied LaTiOa/SrTiOa interface J^JIii^J^ 
For the present interface, the atomic polarization occurs 
mostly in the STO part, with the oxygen atoms moving 
towards the interface, while the cations move away from 
it, so that there is a buckling in the planes. The largest 
cation-anion polarization of about 0.15 A occurs for the 
interfacial Ti02 and the SrO layers as seen from Fig. [51 

The layer-projected densities-of-states (Fig. [3]) calcu- 
lated with the LMTO method shows that the half elec- 
tron per unit cell needed to overcome the polar catas- 
trophe resides on the STO side and spreads out to sev- 
eral Ti layers near the interface. The shift of the band 
edges from layer to layer indicates an electric field dis- 
continuity consistent with the accumulation of a positive 
monopole charge at the interface. Electrons attracted to 
this monopole reside on the STO side, as expected from 
the computed band line-up at the interface (Fig. [3]) and 
confirmed from the charge density contour plot (Fig. [5]). 

An important result of this work is the existence of 
two types of carriers at the interface with possibly differ- 
ent transport properties, which is inferred by examining 
the conduction band structure presented in Fig. [S] The 
occupied conduction bands are made out of Ti(d) t2g 
states. As indicated from the wave function characters 
of Table I, the lowest band is well localized on the first 
Ti layer because of a much lower interface potential there 
(Recall that the polar catastrophe argument would place 
as many as 0.5 electrons per cell at the interface). The 
xy bands are lower than the yz or zx bands because of 
the larger hopping integrals along both directions in the 
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FIG. 3: Layer-projected densities-of-states indicating the 
spread of the interface electrons (0.5 e~/cell occupying Ti(d) 
bands) into several layers near the interface on the STO side. 
The Fermi energy is taken as zero. 




FIG. 4: Computed band offset at the interface. Valence band 
positions were computed by using the 0-ls core state as the 
reference level. The conduction band edges were inferred by 
using the experimental band gaps.— 

xy plane, parallel to the interface. Above these bands 
at the r point are several other bands which occur close 
together in energy. These are xy-type bands located on 
the second or the third Ti layer and the yz or xz-type 
bands, the latter being made out of a linear combination 
of states on several Ti layers and have small dispersions 
along one of the two planar directions. Apart from the 
lowest band, states in the remaining bands spread sev- 
eral layers into the STO part and because of the finite 
size supercell used in our calculations, the bands belong- 
ing to the two interfaces mix, producing a splitting of 
the energies, visible, e.g., near the X point in the figure. 
This particular split becomes amplified in the Fermi sur- 
face, producing the two outermost wings both along X 
and Y directions (inset of Fig. [6|). The two outermost 
wings of the Fermi surface would not split without this 
interaction. 

The basic features of the Fermi surface for the iso- 
lated interface is illustrated in Fig. [71 which is extracted 
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FIG. 5: (Color online) Charge contours of the occupied states in the conduction bands (Fig|6} plotted on the xz-plane (normal 
to the interface). Electrons penetrate very little to the LAO side. Note that only the xz orbitals, but not the xy or yz orbitals, 
show up on such a plane. Nevertheless, the contours faithfully reproduce the spread of some of the electrons (see text) deeper 
into the STO part. The starting density for the contours is 3.2 x 10~^e/A'' and successive values increase by the factor of 3. 
Inset shows the Ti (d) occupancy as a function of the distance of the Ti layers from the interface. 
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FIG. 6: Conduction bands near the Fermi energy and their 
predominant orbital characters. Band dispersions are shown 
along the interface (xy plane) with X — n/a{l, 0, 0) and M — 
7r/a(l, 1, 0), where 'a' is the in-plane lattice constant. 



from the calculated Fermi surface of Fig. [5] and where, 
for pedagogical reasons, we have omitted the hybridiza- 
tion between the bands. The circles are derived from the 
xy orbitals, which disperse with the tight-binding ma- 
trix element of VddTT along both x and y, resulting in the 
tight-binding band structure: €^^'^\k) = 2Vdd7T {cos kxa + 
cos ky a), which is isotropic for small Bloch momenta, 
while the xz and the yz orbitals disperse with the tight- 



TABLE I: Square of the wave function characters for typ- 
ical conduction states at the F point, labelled in Fig. 
|6l The Ti02 planes have been numbered consecutively 
starting from the interface. 



State (Ti02)i 


(Ti02)2 


(Ti02)3 


(Ti02)4 


rest 


i>i 0.82 


0.02 








0.16 


i/)2 0.02 


0.17 


0.44 


0.26 


0.11 


ips 0.01 


0.35 


0.18 


0.37 


0.09 
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Til(xy) ; ; '^^Til, Ti2, Ti3;xz) 



Ti1,Ti2, Ti3;yz) 



7t/a 



Ti2, Ti3, Ti4(xy) 



FIG. 7: (Color online) Idealized Fermi surface for the iso- 
lated interface, extracted from the DFT results of Fig. |5] 
Hybridization between bands rounds out the sharp crossing 
points in the Fermi surface, leading to more complex orbits 
as seen in the inset of Fig. [S] 



binding matrix element VddTr along one direction but with 
a much smaller^^ VddS along the other, so that we have 
for the xz bands: e(^^)(fc) ~ 2VddTT cos kxa+2VddS cos kya, 
and similarly for e(^^-'(fc), which produce the elliptic 
Fermi surfaces. 

Our results suggest a possible explanation for the ob- 
served low carrier density for the samples grown under 
oxygen rich conditions. This growth condition is thought 
to produce an intrinsic interface, where the observed car- 
rier density is ~ 1 — 2 x 10^^ cm^^, corresponding to only 
about 0.015-0.03 electrons per interface unit cell, much 
smaller than the half electrons per cell needed to sup- 
press the polar catastrophe. These electrons need not 
reside on just one interfacial layer, rather, they can be 
spread over several layers and the polar catastrophe is 
avoided so long as the total electron count is right. Not 
all of these electrons might contribute to transport. 

From Fig. [7| we infer the presence of two types of 
electrons at the interface. The first originates from the 
lowest conduction subband, which has a strong 2D char- 
acter consisting mostly of Ti(d) states on the first layer 
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and forming the largest circle in the Femi surface plot. It 
is well-known that all states in 2D are Anderson localized 
by disorder and these states may therefore become local- 
ized. In addition, these electrons may form self-trapped 
polarons due to the Jahn- Teller coupling of the Ti(d^) 
center. The unpaired electron could then behave as a 
Kondo center. The xz/yz states which form the ellip- 
tic Fermi circles in Fig. [7] are also prone to localization 
on account of their heavy masses along a planar direc- 
tion. All these electrons may therefore not participate in 
transport. 

The second type of electrons are from the Ti (xy) bands 
above the lowest band which spread over several Ti lay- 
ers into the bulk and show up as the three innermost 
circles in the Fermi surface plot. These electrons have 
small masses along the plane, they spread several Ti lay- 
ers into the bulk and are as a result difficult to localize by 
disorder or electron-phonon coupling and we argue that 
they are likely the mobile electrons that would contribute 
to transport in the intrinsic samples. Since the electron 
density in the Ti(d) bands is low, far from half filling, 
electron correlation effects are not expected to alter the 
metallic behavior apart from a mass renormalization in 
a Fermi liquid picture. 

The number of these mobile electrons, estimated from 
the area of the three innermost circles in Fig. [7] is about 
8 X 10^'^ cm^^, which is five times smaller than the elec- 
trons needed to satisfy the polar catastrophe. It is how- 
ever still several times higher than the measured car- 
rier density in the intrinsic samples. This discrepancy 
could be due to several reasons. Density-functional cal- 
culations tend to yield too high an energy for localized 
states, which would mean that the density in the delocal- 



ized bands is overestimated. Another possibility is that 
the mobile carrier density may be reduced due to addi- 
tional localization effects not considered here explicitly. 
So, according to this scenario, the mobile electrons arise 
from the inner pockets in the Fermi surface, while the 
remaining electrons are localized. 

There are indeed signatures of these two types of car- 
riers in the experiments. For instance, Brinkman et ali^ 
were able to explain their transport data in terms of two 
types of carriers. These experiments showed a Kondo 
resistance minimum as a function of temperature, which 
they interpreted in terms of itinerant electrons scattering 
off of localized Ti (3d) moments. Within our scenario, the 
formation of the local moments is entirely likely due to 
the localization of some of the electrons. The delocalized 
electrons lead to the possibility of superconductivityi^ 

It would be valuable to measure the Hall effect on the 
intrinsic samples, where several authors have been un- 
successful in observing the Shubnikov-de Haas (SdH) os- 
cillations. It is possible that the oscillations become sup- 
pressed due to the intermixing of the localized and the 
itinerant parts of the Fermi surface, leading to strong 
scattering between them before the completion of full 
magnetic orbits. In that case, one might have to go 
to higher magnetic fields, where magnetic tunneling be- 
tween the bands might lead to complete orbits. 
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